Bulk Inversion Asymmetry effects on the band structure of zincblende 
heterostructures in an 8-band Effective Mass Approximation model 
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We have developed an 8-band Effective Mass Approximation model that describes the zero field 
spin splitting in the band structure of zincblende heterostructures due to bulk inversion asymmetry 
(BIA). We have verified that our finite difference Hamiltonian transforms in almost all situations 
according to the true D2d or C2v symmetry of [001] heterostructures. This makes it a compu- 
tationally efficient tool for the accurate description of the band structure of heterostructures for 
spintronics. We first compute the band structure for an AlSb/GaSb/AlSb quantum well (QW), 
which presents only BIA, and delineate its effects. We then use our model to find the band struc- 
ture of an AlSb/InAs/GaSb/AlSb QW and the relative contribution of structural and bulk inversion 
asymmetry to the spin splitting. We clarify statements about the importance of these contributions 
and conclude that, even for our small gap QW, BIA needs to be taken into account for the proper 
description of the band structure. 

PACS numbers: 73.21. Cd,73.21.Fg 
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I. INTRODUCTION 



In recent years, integeat in developing spin-sensitive 
devices (spintronics)El'E(Era has fueled renewed investiga- 
tions into spin phenomena in semiconductors. The aim 
is to control not only the spatial degrees of freedom of 
the electron, but also the spin degree of freedom. Use- 
ful spintronic devices can be devised if such control is 
achieved j-||A,Humber of such devices have already been 
proposedBfl'QH. If a full understanding of the operation 
of spintronic devices is desired, a thorough knowledge 
of the band structure including all spin details will be 
needed. In particular, detailed knowledge of heterostruc- 
ture bands is required when studying, and perhaps tai- 
loring, time evolution and transport phenomena of spin 
ensembles. 

Most of the standard implementations of the Effective 
Mass Approximation (EMA) for III-V heterostructures. 
do not take into account the bulk inversion asymmetrytl 
(BIA) present in zincblendesOliil. The BIA lifts Kramers 
degeneracy and, therefore, a potentially important source 
of spin splitting in the bands is not included. There have 
been several proposed modifications to EMA to account 
for BIjL^nd its associated spin splittings, ranging from 
2-bandy to-J4-bandy and IG-bandllJ models. Also, Zhu 
and Changll3 started from an 8-band model to generate 



perturbatively a 2-band Hamiltonian for electrons and a 
4-band Hamiltonian for holes, and they performed their 
calculations of inversion asymmetry eflEects in that re- 
duced basis set. Rossler, Winkler et al&jtl have used an 
8-band model solved with a quadrature method that in- 
cluded BIA to study Al^jGai-ajAs/GaAs quantum wells, 
where the spin splitting is much smaller than the one 
shown here. In this work, we study the symmetry prop- 
erties of the heterostructure Hamiltonian and verify that 
the 8-bajid model can account for both the BIA and 
Rashballa effects on the spin splitting. 

We present a Finite Difference Method (EDM) imple- 
mentation of a compact 8-band EMA method capable of 
reproducing the D2d and C2v symmetry of [001] grown 
quantum wells and superlattices. We have used this 
method to study the relative contribution of structural 
inversion asymmetry (SIA) vs. BIA on the spin split- 
ting in the conduction subbands. Our study has allowed 
us to put into perspective the, apparently contradictory 
statements of Lommer et altB, stating that BIA (SIA) 
dominates the splitiirig in large (narrow) gap materials, 
and Cardona et a/.cfl, listing a proportionality constant 
for the splittingEI that becomes roughly larger with 
smaller bandgap (see Table |). The EDM allows for ar- 
bitrary external fields and the dcscijiotion of tunneling 
phenomena with only a few changesEj, while maintain- 
ing numerical efficiency. 
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Section || describes how the 8-band bulk Hamiltonian 
yields the EMA equations and the method for their so- 
lution. In Sees. [II and IV the effects of bulk inversion 



has been defined as 



asymmetry in symmetric and asymmetric quantum wells 
are explored, taking as example the band structure of 
AlSb/GaSb/AlSb and AlSb/InAs/GaSb/AlSb quantum 
wells. Finally, the results are summarized in Sec. 



II. EMA HAMILTONIAN WITH BIA 

For the calculation of the heterostructure bands, the 
Effective Mass Approximation (EMA)E3 based on an 8- 
band k • p formalism is used, Jlbere are several published 
8-band k • p Hamiltoniana23'EilO, each including a more 
or less detailed set of effects. Here_we take the Hamilto- 
nian constructed by Trebin et alrB as a starting point. 
This Hamiltpniap is constructed by means of an invariant 
expansioncJ'LZlH and, when applied correctly, guarantees 
the inclusion of all matrix elements compatible with the 
Td symmetry group of the zincblendes up to the desired 
order in the electron wavevcctor k. 

For completeness. Appendix |^ reproduces the full 8- 
band k • p Hamiltonian from Ref. ^ Note that, due to 
the way that it has been constructed, this Hamiltonian 
takes into account all the effects of the spin-orbit inter- 
action in the matrix elements up to k^, and in particular 
the second order s—p coupling via remote Fs statesEj re- 
sponsible for most of the contribution to the spin splitting 
in the conduction band. Strain and coupled strain/spin- 
orbit effects are also properly described by this method. 
Note that Appendix |A| corrects two typographical errors 
present in Ref. |2^. 

The 8-band Luttinger parameters ji apeearing in Ta- 
ble ^ are the true Luttinger parametersE^I with the 
conduction band contribution subtracted becawse its ef- 
fects are treated exactly. They are related byCJ 
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where Eg is the energy gap of the compound, and Ep 



Material 


Band Gap (eV) 


7c (eV-A^) 


GaAiS 


a 


1.52 


25.5 


InPt 




1.42 


8.5 


GaSb 




0.81 


186.3 


InAs 


b 


0.418 


130 


InSbl 


] 


0.235 


226.8 



"Adapted from Ref. F3. 
'Ref.^ ^ 



Et 



(2) 



with TTT-e being the free electron mass and P the irre- 
ducible momentum matrix element. 

Table H shows the numerical values used in the cal- 
culations for the parameters of a number of materials. 
B and C are parameters describing the BIA effects. As 
shown below, B is mainly related to the conduction band 

splitting while C is the coefficient for the linear k split- 
ting in the valence band. No magnetic field effects will 
be included. For structures grown along the [001] direc- 
tion, C2, C4 and C5 will not be needed because there 
is no shear stress. Finally, the effect of remote bands 
on the conduction effective mass will be neglected. This 
amounts to setting the A! parameter to zero. 



Analytical expressions for the bands close to 
the zone center 



Starting from the full Hamiltonian shown in Ap- 
pendix Q analytical expressions for the bands near the 
zone center can be found. These are useful for finding 
measurable quantities such as effective masses and intra- 
band splittings as a function of the model parameters. 
These expressions can also be useful when relating the 
parameters of the model to the parameters used in other 
families of k • p Hamiltonians. The effective masses as a 
function of the parameters can be found, for example, in 
Ref. ||. 

Along [100], the dispersion relation to lowest order in 
k for the heavy hole (HH) band is 



Ehh {\^)=E,+Ck^, 
while for the light hole (LH) band 

Elh (k) =E,- Ck,, 



(3) 



(4) 



TABLE I: Band Gap and 7c for selected III- Vs. 



where E^ is the energy at the zone center. Note that each 
band along [100] is still doubly degenerate as required by 
the Td symmetry. In this case C describes a linear split- 
ting between the heavy and light holes. As a consequence 
of BIA the top of the valence band lies slightly .away from 
fc = 0. This has been experimentally observedujO. 

Figure |l| shows the effect of the inclusion of the param- 
eter C in the HH and LH bands along [1 00] calculated 
from the numerical diagonalization of Eq. (Al). Clearly, 
the effect is not appreciable at normal scales of k and 
energy. As shown in the inset, we need to look very close 
to the zone center to observe the effect of C. Note the 
very small value of the energies involved. Away from the 
zone center the bands with BIA effects recover the usual 
behavior except for some additional splitting present. 

Along the [110] direction, the energy of the electrons in 
the CB as a function of the wavevector k is, up to third 
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Parameter 


InSb 


GaSb 


AlSb 


InAs 


GaAs 


AlAs 


a {ky 


6.4794 


6.096 


6.136 


6.058 


5.653 


5.66 


E, {eVf 


0.235 


0.813 


2.219 


0.356 


1.52 


3.002 


{eVYr 





0.56 


0.11 








-0.55 


^so 


fV)l 


0.803 


0.8 


0.75 


0.41 


0.341 


0.279 


1 




2.59 


2.58 


1.44 


2.05 


2.01 


1.74 


T 




-0.6 


-0.58 


-0.35 


-0.44 


-0.41 


-0.37 


1 




0.67 


0.65 


0.39 


0.48 


0.46 


0.42 


C {eV-k)^ 


-9.32x10"^ 


7.00x10"'' 





-1.12x10"^ 


-3.40x10"^ 


2.00x10 


P {eV-k)% 


9.35 


9.21 


8.41 


9.17 


9.86 


8.94 


B [eV-k^Y 10.3 


49.9 





13.7^ 


30.4 


21.3" 


Ci [eVf 


-6.17 


-6.85 


-6.97 


-5.08 


-7.17 


-5.64 


Cii (GPa) 


69.18 


88.34 


87.69 


83.29 


112.6 


120.2 


Ci2 (GPa) 


37.88 


40.23 


43.41 


45.26 


57.1 


57 


C44 (GPa) 


31.32 


43.22 


40.76 


39.59 


60 


58.9 



"Rcf. 21 



"Rcf. 
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ence band offsets are consistent within the systems comprised of (InSb), (GaSb, AlSb, InAs) and (Al^jGai-^! As). 



"^The va 
''Ref. |l^ 

■^From 7c obtained in Ref. 23 
■fProm 7c obtained in Ref. 21 
fProm 7c obtained in Ref. T2 

''Adapted from the Bir-Pikus deformation potentials a, b, d in Ref. 



TABLE II: Parameter values for some materials. 
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k II [100] [10"" A'] 



-0.05 0.00 0.05 

Ml [100] [A-^] 

FIG. 1: HH and LH bands along [100] near the zone center 
with (BIA) effects included. No difference in the behavior 
of the bands is appreciable at this scale. The inset shows a 
blowup of the region very close to the zone center. It is seen 
that the inclusion of C 7^ induces the presence of a linear 
splitting in the vicinity of k = 0. Note the scale of the inset. 



order, 

£;cB(k) = i^, + -J^±i7cfc^ (5) 

with ra^jQ being the conduction band effective mass, and 
the I? sphtting coefficients 7c given, in terms of the 



model parameters, by 

_ P 2BEgAso - V3CP {Eg + Agp) 
' 3 (Eg + Aso) 

It is easy to shovif that the contribution to 7c of the 
part containing C is only about 4% for InSb and InAs. 
That contribution goes down to about 0.3% for GaAs and 
AlAs, and it drops to a mere 0.03% for GaSb. Therefore, 
it is a good approximation to consider that all the split- 
ting in the conduction band is due to the nonvanishing 
bulk inversion parametec-jB. Then, if C is neglected, the 
result of Eppenga et alEj for the expression of 7c is re- 
covered. Note that, in order to turn off BIA effects, both 
parameters B and C need to be set to zero. 

As in the [100] case, the inclusion of C changes the 
characteristics of the bands very close to the F point. 
Again, it provides them with a small linear component. 
But, in contrast to the [100] case, here the HH and LH 
bands are not doubly degenerate. The linear splittings 
^LH [110] and Affff ^hq^ turn out to be the same for both 
HH and LH bands, and are given by 

^HH [110] = ^LH [110] = ViCk. (7) 

This result is slightly different from the one indicated 
in Eq. (7.5) of Ref. Q A numerical diagonalization of 
the Hamiltonian has been performed to check the validity 
of Eq. ( [^) The discrepancy arises because the splittings 
in Ref. ]ij are valid in the region where the quadratic 
(effective mass) splitting predominates, while the result 
obtained here is valid in the region where the linear split- 
ting dominates. 
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FIG. 2: Band structure for GaSb near F along [110] and spin 
splitting of the conduction band. The dashed line of plot b) 
corresponds to the function Splitting — "/ck'^ , with 7c = 186.3 
eV-k^ 



The SO band also presents splitting A^q [ho] 7 pro- 
portional to the B parameter only: 

2BP 



A 



so [110] 



3iEg + Aso) 



(8) 



Along [111], in the region where the linear splitting 
dominates, the heavy hole (HH) band has the dispersion 
relation 



Ehh Ck)^E,±V2Ck~ 



2l2 



2m, 



(7iL- 273L)+0(fc*), (9) 



while for the light holes 

Elh (k) = E„-^ (7ii + 273l) + 0{k^). (10) 

The light hole, conduction and split-off bands are de- 
generate along the [111] direction, as can also be deduced 
by group theory argumentaj. The linear splitting ob- 
tained for the heavy holes agrees with the result of Car- 
dona et aZo. 



B. Bulk bands and spin behavior 
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FIG. 3: Direction for the spin of the spin-split states of the 
lowest conduction subband of GaSb. This plots sweeps a 
circular path in k space with fcz = 0. No spin direction is 
specified for the (100) family because the states are spin de- 
generate. 



combinations of components of k and constructing an 
invariant 2-band Hamiltonian for the conduction band, it 
is found that the Hamiltonian wilLinclude the following 
term breaking the spin degeneracjciJ: 

Hki = Ic [o'xkx {ky — /cf) -I- cyclic permutations] . (11) 

The two band Hamiltonian in Eq. ([l^) also predicts 
the direction where the split spins will be pointing. For 
example, it is easily seen from the previous equation that 
if k = (1,(5,0), with 8 a positive infinitesimal, the spin 
will point along the ±?/ direction. Similarly, symmetry 
requires the spin to point along (—1, 1, 0) or (1, —1, 0) for 
k along the [110] direction. This is indeed obtained in 
the numerical diagonalization of the Hamiltonian includ- 
ing the BIA effects, as seen in Fig. |[ That figure shows 
a circular sweep in k space with fc^ = 0. The arrows rep- 
resent the direction towards which the spin of the lowest 
conduction subband states is pointing. The horizontal 
axis represents k^^ while the vertical axis indicates the 
ky component of the state. The states belonging to the 
(100) directions are spin degenerate; therefore, no spin 
direction is given for them in Fig. |^. 



An example of bulk bands computed using the Hamil- 
tonian with the full symmetry can be seen in Fig. Ha), 
which shows the band structure of GaSb along the [110] 
direction. For this direction, all bands are spin split ex- 
cept at the F point. Plot b) shows the energy splitting of 
the CB states as a function of fc. Close to the zone cen- 
ter, the splitting follows the behavior described in Eq. (||), 
with 7c taken to be the experimental value. In this case, 
it is seen that the usual 2-band Hamiltonian model [see 
Eq. (|ll]) below] can describe the splitting with good 
accuracy up to 2.5.% of the zone boundary. 

This splitting^ can also be predicted by the methods 
described in Refs. |2^,^. Going up to order 3 in the 



Interface conditions and Hermiticity in the 
EDM 



In Appendix ^ we describe the Finite Difference 
Method (FDM) we use to find the heterostructure bands 
and eigenstates. This method transforms the EMA set 
of coupled ordinary differenti al e quations into the eigen- 
value problem shown in Eq. (C9). 

The hermiticity of the discretized Hamiltonian oper- 
ator in Eq. 



dCQl) is ensured if H^_^i = Hi+i,,. bmce 
the H'--')'s are themselves Hermitian, an inspection of 
Eqs. ( |C10[ )-( |CT^ ) shows that this is indeed the case. It 
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From this it must be concluded that the enforcement of 
interface conditions of the form ( |l^ ) is not possible if the 
Hermiticit y of t he Ha miltonian is to be preserved, and 
that Eqs. ( |C1C| )-( C12 ) must be used instead for all mesh 
points. 



12 3 4 

FIG. 4: Mesh used in the study of interface conditions. 



is also clearly seen that the introduction of quantum well 
or superlattice boundary conditions (BCs) doesn't affect 
the Hermiticity of the Hamiltonian. 

There exist in the literature several proposals on what 
are the correct qua atitigs . to . match at the interface be- 
tween two materialsEaooO. Most of them require the 
continuity of the envelope function and a quantity that 
has the general form: 



[Ad, + B] F 



(12) 



where the matrices A and B take different values depend- 
ing o n th e aut hor. Using the finite difference formulae 
Eqs. (C6)-(|C7|), the conti nuit y of Eq. ( |l^ ) can be written 
in a form similar to Eq. (C8): 



H2,lFi + H2,2F2 + H2,3F3 







(13) 



where the subindexes i, j are referred to the mesh points 
in Fig. ^ and the H^j are given by 



H2,l = +1- 
H2,2 = B^ 

H2,3 = 



4Az 
B« 

- A« 



4A;; 



(14) 



with L (R) meaning the material at the left (right) of the 
interface. Isolating F2 from Eq. (|l^) and plugging it into 
the corresponding equations for the F^'s, one obtains 



Hi^oFo + (Hi^i 

Hl,2H2 2H2,3F3 

H 



-H3_2H2 2H2,lFi 



H3.2H ^ 



Fi- 



lls) 



3,3 — -n-3, 2-1-12,2 
H3_4F4 — EF4 



^2,3) F3 + 



(16) 



Now, in order to preserve the Hermiticity of the dis- 
cretized Hamiltonian, one should have 



H 



t 

2.3 



H 



3,2 



H 



t 

2,1 



Hi 



(17) 



However, looking at Eqs. (0) and E^s. (|cTo|)-(|cT2|) 
one sees that the requirements in Eq. ( |l7[ ) are not satis- 
fied by the discretized version of the interface conditions 



D. Symmetry group of the discretized Hamiltonian 

For the rest of this article, we will refer to "symmet- 
ric" and "asymmetric" heterostructures, mainly quan- 
tum wells (QWs). Here, "symmetric" will be taken to 
mean that the sequence of materials and their respective 
thicknesses are left unchanged under the inversion oper- 
ation (i.e., they are macroscopically symmetric). Thus, 
an AlSb/InAs/AlSb QW is called symmetric, while an 
AlSb/InAs/GaSb/AlSb QW is called asymmetric. This 
definition is made in order to avoid confusion with the mi- 
croscopic symmetry -that is, the symmetry group- of the 
QW. All asymmetric [001] heterostructures made from 
zincblendes are described by the point group. On 
the other hand, symmetric [001] QWs can belong to ei- 
ther the C2V or the I?2d symmetry groups depending on 
an interplay of characteristics such as the parity of the 
number of monolavejrs, the existence of a common atom 
in the constituenttcS or the presence of a supplementary 
half-layercj. 

T he sy mmetry group of the discretized Hamiltonian in 
Eq. ( |C9| ) can be found by verifying that all operations g 
of the Td point group satisfy the relationship^ 



D-\g)Ji(kn)Dig)=llig-\ 



(18) 



where D(g) is the representation of the g operator in 
the basis of the discretized Hamiltonian H(k||), and 
k|| = {kx, ky). This tedious procedure can be done wilh 
the help of computer software, such as MathematiccsM, 
which automates algebraic manipulations. It is seen that 
the EMA Hamiltonian corresponding to symmetric struc- 
tures transforms according to D2d, while for asymmetric 
structures it transforms according to C2V ■ This is in con- 
trast to the majority of EMA implementations, which 
lack the inclusion of bulk inversion asymmetry effects and 
reproduce an approximate D411 symmetryEJ for symmet- 
ric structures. 

Tables HI and |^ show the requirements that the un- 
derlying symmetry of the atom arrangement imposes on 



Point/Line 


Point Group 
Symmetry 


Spin 
Splitting 


Spin directions 


r 




No 




A [100] 


C2 


Yes 


[100], [100] 


E [110] 


Cs 


Yes 


[110], [110] 


Other points 


Ci 


Yes 


Undetermined 



TABLE III: Symmetry requirements on spin splitting and di- 
rections for points in the — ky plane in a D2d structure. 
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FIG. 5: Layer arrangements for a no-common-atom quantum 
well. An arrangement as in a) in the EMA would yield D2d 
symmetry. The alternative arrangement b) yields the correct 
C2v symmetry of the heterostructure. 



at the interface are the same, there is an asymmetry in 
the bond orientation, which the present EMA method 
cannot take into account. 

With the control that the proposed model allows over 
the symmetry of the Hamiltonian, it is possible to use 
the EMA in studies of the spin splitting appearing in 
heterostructures. Thus one can delineate the role of 
bulk inversion asymmetry vs. structure inversion asym- 
metry, layer asymmetry vs. interface asymmetry, etc. 
This model also provides a straightforward and easy to 
implement tool for studying reduced sym pietr y effects, 
such as the presence of optical anisotropymEJ, and the 
mixijag.of heavy hole and light hole states at the zone cen- 
teiliOo in QWs and superlattices. Some of these works 
have followed the alternative approach of introducing in- 
terface parameters to describe the lower ing of symmetry 
respect to the standard EMA formalismc§E3E3. In these 
cases interface asymmetry effects are well described, but 
BIA is not accounted for. 



the spin degeneracy of the energy levels and the direc- 
tion where the spins are pointing in case the levels are 
not degenerate. Note that, for a structure with D2v sym- 
metry (|i-f., BIA effects only) the bands along [100] must 
be splitli3 even though there is no splitting along [100] 
for bulk zincblendes. 

Although the k • p method is not designed taking into 
account the interface characteristics of no-common-atom 
(NCA) heterostructures, it is possible in some situations 
to modify the simulated structure to obtain at least the 
right symmetry effects. Figure |^.a) shows a NCA quan- 
tum well. The way that the boundaries of the layers are 
set up, the well would be symmetric and, therefore, the 
Hamiltonian would have D2d symmetry instead of the 
C2v corresponding to the asymmetric interface bonds. 
However, the material boundaries in the k • p method 
are arbitrary to half a monolayer. As seen in Fig. ^.b), 
a simple rearrangement of the material boundaries re- 
produces the asymmetry in the bonds and allows us to 
take into account, at least qualitatively, the effects of the 
lower symmetry. 

The only case that cannot be modeled through these 
rearrangements is when a common atom structure, such 
as an AlAs/GaAs/AlAs QW, has C2v symmetry due to 
the well having an odd number of monolayers. In that 
case, even though the species participating at the bond 



Point/Line 


Point Group 
Symmetry 


Spin 
Splitting 


Spin directions 


r 




No 




A [100] 


Ci 


Yes 


Undetermined 


E [110] 


Cs 


Yes 


[110], [110] 


Other points 


Ci 


Yes 


Undetermined 



TABLE IV: As in Table III, for a C2v structure 



III. BULK INVERSION ASYMMETRY 
EFFECTS ON SYMMETRIC QUANTUM WELLS 

In this section, the methods described in Sec. || will 
be used to calculate the electronic properties of a sym- 
metric quantum well. In particular, focus will fall on 
AlSb/GaSb/AlSb quantum wells. However, some of the 
results derived are a consequence of the underlying sym- 
metry of the structure rather than the constituents them- 
selves. Therefore, these particular results will illustrate 
general considerations. 



A. SQWs without BIA terms 

Figure ^ shows the band structures along the [100] and 
the [110] directions of a common atom AlSb/GaSb/AlSb 
symmetric quantum well (SQW) grown along the [001] 
direction and with a well thickness of 8 monolayers (24.4 
A). Since no inversion asymmetry affects are included, 
the bands show Kramers degeneracy throughout the Bril- 
louin zone and the quantization axes of the spins are not 
unambiguously defined. 

The labels El, HHl, LHl and HH2 shown in the plots 
correspond to the first electron, first heavy hole, first light 
hole and second heavy hole states in the QW respectively. 
They refer to the main bulk state contribution at k = 0. 
For a well without BIA terms and in the zone center, 
the heavy holes decouple from the rest of the bands, and 
the HHn states have only bulk heavy hole components. 
This is in contrast to the En (LHn) bands, which have 
small bulk light hole (electron) and split off contributions 
even at the zone center due to the loss of translational 
symmetry along [001] caused by the well potential. 
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FIG. 8: a) Spin splitting dependence for an AISb/GaSb/AlSb 
QW along the [100] line. The calculated proportionality con- 
stant is asiA = 22 X 10"^'' eV-cm. b) Splitting dependence 
along a circle in the — ky plane, with k = 0.001 A. 



FIG. 6: Bands along [100] and [110] for an AlSb/GaSb/AlSb 
SQW 8 monolayers thiclt without BIA terms. 
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FIG. 7: Same as Fig. E but with BIA terms. 



B. SQWs with BIA terms 

Figure ^ shows the same band structures as in Fig. ||, 
but with BIA terms included. As predicted by group 
theory (cf. 



Table [II 



the bands are split along both 
directions except at the zone center. This is the major 
difference with most of the EM A models in the literature. 

With the inclusion of BIA terms, the heavy hole states 
couple with the light holes by means of remote states 
through a perturbative mised spin orbit and k • p inter- 
action parametrized by Co. Thus, the HHn states lose 
their pure bulk heavy hole character and, in particular, 
the HHewen (HHorfd) mix with the LHodd (LIIe?;en)E3. 
However, looking at the wavefunction for the HHl state 
of the structure in Fig. |^, it can be seen that the contri- 



bution from bulk components other than the HH to the 
probability density is about 8 orders of magnitude less 
than the heavy hole contribution. 

Figure ^ shows the linear behavior and the isotropy of 
the conduction subband spin splitting close to the zone 
center. Plot a) shows the dependence of the spin split- 
ting in the CB along the [100] line. It is seen that the 
splitting is linear close to the F point, with a proportion- 
ality coefficient of asiA = 22 x 10^^" cV-cm, where the 
splitting is 



^BIA = 2aBIAk. 



(19) 



Although there is an isotropic linear spin splitting, this 
must not be confused with the Rashba splittinglj. By- 
chkov and Rashba introduced a splitting coefhcient an 
in the context of asymmetric quantum wells. The split- 
ting studied there is derived from a model Hamiltonian 
that describes only structural inversion asymmetry (SIA) 
effects, but not bulk inversion asymmetry. As a conse- 
quence, the spin directions that they predict don't apply 
to the SQW situation (cf. Sec. 0). 

The computed aBiA for this structure is about half 
of some of the highest-^jcedicted Rashba coefficients for 
asymmetric structuresc3E2l where only SIA contributions 
are taken into account. This shows that BIA effects need 
to be carefully studied before neglecting them in a calcu- 
lation. 

A model Hamiltonian for spins in the conduction sub- 
bands |Of SQWs in the same spirit as the RaaUpa Hamil- 
tonianta was obtained by de Andrada e SilvaEl 



H = as I A {cTxkx 



CTyky) 



(20) 



where a are the Pauli matrices, olbia — 
- {Fs;kx,ky\^cdl\Fs;kx,ky) and \Fs;kx,ky) is the 
envelope function corresponding to the electron traveling 
in the plane with wavevector kx,ky. 

The behavior of the spins when BIA terms are included 
is very interesting. Figure |^ shows the spin directions of 
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FIG. 9: Spin directions for the lower conduction subband of 
an AlSb/GaSb/AlSb SQW. The thickness of the weU is 8 ML 
(24.4 A). The spins are plot at 15 ° intervals, and correspond 
to states lying on a circle in the kx — ky plane with k = 0.01 

A. 



the eigenstates of the lowest conduction subband from 
our calculations. The spin directions are shown for a cir- 
cular sweep in the — ky plane keeping k = 0.01 A. The 
directions of the spins agree with what would be pre- 
dicted from Eq. (^0|). The spins at a given point in the 
kx — ky plane point in opposite directions for the two sub- 
bands. Note that in a given subband, although the x and 
y axes are equivalent, in one of the axes the spin points 
outward while in the other it points inward. The expla- 
nation lies in the way that the x and y axes are connected 
and in the fact that spinors don't change sign under in- 
version. For a QW with D2d symmetry, the x and y axes 
are equivalent through a reflection by the plane contain- 
ing the [110] and [001] directions. The reflection by this 
plane can be thought of as a rotation of 180° along the 
[110] direction followed by an inversion. Starting with a 
state with ky component only \ky, ty) (spin pointing out- 
ward), the rotation will send it to | — kx, T-i) (still out- 
ward). Then, the inversion will flip k, but not the spin, 
sending the state to \kx, T-s) (spin pointing inward). 



IV. BULK INVERSION ASYMMETRY 
EFFECTS ON ASYMMETRIC QUANTUM 
WELLS 



In this section the structure under study will be 
an AlSb/InAs/GaSb/AlSb asymmetric quantum well 
(AQW) grown along the [001] direction compliant with 
a GaSb substrate. The thickness of the InAs and GaSb 
layers is 8 monolayers (ML) each, with a monolayer thick- 
ness of 3.048 A. 
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FIG. 10: Bands without BIA effects for an 
AlSb/InAs/GaSb/AlSb AQW grown along the [001] di- 
rection compliant with a GaSb substrate. The bands are 
along [100] and [110]. The thickness of the InAs and GaSb 
layers is 8 ML each. 

A. AQWs without BIA terms 

Figure shows the band structure along the [100] and 
the [110] directions of the AQW without the inclusion of 
BIA terms. The structural inversion asymmetry (SIA) 
reduces the symmetry group from D^h to C2v , and a spin 
splitting appears between the conduction subbands. The 
splitting due to SIA effects is usually modeled using-a 
Hamiltonian first introduced by Bychkov and Rashbali3: 

Hr = ttfl (cr X k) • iv, (21) 

where an is the Rashba coefficient, cr is a vector com- 
posed of the Pauli matrices and v is the axis of symmetry 
of the structure. This Hamiltonian is valid for describing 
the SIA contributions close to the zone center. It predicts 
a linear and isotropic splitting 

Ar = 2aRk, (22) 

where k is the magnitude of the electron wavevector. It 
also predicts that the spins will point tangentially to the 
circles of constant k in the kx — ky plane, which is verified 
by numerical calculations (see Fig. [III ). 

B. AQWs with BIA terms 

The band structure for the asymmetric quantum well 
(AQW) under study with BIA effects is shown in Fig. |lj. 
The effects of inversion asymmetry are highly anisotropic 
in bulktl, and this reflects on the directional dependence 
of the bands. Comparing with Fig. |o|, it is seen that 
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FIG. 11: Spin directions for the lower conduction subband 
of an AlSb/InAs/GaSb/AlSb 8/8 AQW without BIA effects 
included. The spins are plot at 15 ° intervals, and correspond 
to the lowest conduction subband states lying on a circle in 
the kx — ky plane with k — 0.01 A. 
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FIG. 12: Bands along [100] and [110] for an 
AlSb/InAs/GaSb/AlSb AQW with BIA effects. 



the BIA effects are necessary to obtain accurate bands 
in the [110] direction. In particular, note that the LHl 
spHtting along [110] is much smaller when BIA effects are 
included. However, we should point out that the LHl 
splitting, with BIA, along [110] becomes much larger due 
to the inequivalency of the [110] and [110] under the C2v 
point group. 

The inclusion of the Hamiltonian ( po| ) keeps the analy- 
sis of the combined BIA and SIA effects quite simple. For 
a [GDI] structure, the SIA and BIA contributions to the 
splitting can be described by a Hamiltonian Hia made 



90 




270 



FIG. 13: Angular dependence of the spin splitting for an 
AlSb/InAs/GaSb/AlSb AQW. The solid line is the 8-band 
model numerical result. The dashed line is a fit using Eq. plj ) 
with OR = 40.3x10"^" eV-cmand obia = IS.OxlO"^'^ eV-cm. 

from the addition of Eq. @) and Eq. (|l|): 

HiA = obia (cTxkx - (Jyky) + an (a^ky - ayk^) = 

Ox {anky + OL^ihkx) ~ Oy (aeiAfcy + OLBk^) , (23) 

where aeiA (afl) is the coefficient describing BIA (SIA) 
effects. From here, making an analogy with the Zeeman 
splitting, it is easy to find that the splitting in the con- 
duction band (CB) close to the zone center will be 



AiA ^2k^a\ + 2aRaBiA sin 20 + , (24) 

where 9 is the in-plane polar angle. 

A full 8-band numerical calculation of the splitting 
along a circle in the kx — ky plane and the 2-band pre- 
diction from expression (|2^) are shown in Fig. |l^. The 
values from the analytic expression show very good agree- 
ment with the numerical results from our 8-band imple- 
mentation. The numerical results are fitted with asiA — 
40.3 X IQ-io eV-cm and ckbia = 15.0 x lO"!" eV-cm. 
This way the BIA effects are quantified, and it must be 
concluded that they need to be taken into account for 
an accurate description of the bands. This is clearly so 
in the [110] direction, where the contributions are added 
linearly, but it is also true in a lesser degree in the [100] di- 
rection, where the contributions are added quadratically. 
This also has a clear effect on the theory of extraction of 
the Rashba coefficient from Shubnikov-de Haas measure- 
ments in two-dimensional electron gases. (2n.F.Gs) , where 
BIA effects have been usually neglectcdLnHEj. 

For a quantum well where the 2-band model is valid, 
the BIA splitting coefficient for the CB can be estimated 

witi£3 

"BIA ~ Y2~' \^^) 

where "fcw and Liy^ are the k^ splitting coefficient of the 
CB and the thickness of the layer where the electrons are 
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FIG. 14: Spin directions for the lower conduction subband 
of an AlSb/InAs/GaSb/AlSb AQW with BIA effects. Tfie 
spins are plot at 15 ° intervals, and correspond to the lowest 
conduction subband states lying on a circle in the — ky 
plane with k = 0.01 A. 



confined respectively. This estimate will become more ac- 
curate as the well becomes thicker. From this expression 
and Eq. (^) it is readily seen that BIA effects will be con- 
siderable when the material in the well layer has a low 
bandgap and high spin-orbit interaction, such as InAs, 
GaSb and InSb (cf. Table ^ and in a narrow well. So, it 
has been seen that, a priori, it is not possible to consider 
only SIA effects for an asymmetric structure even if the 
constituents are kiw bandgap materials. The statement 
of Lommer et al£3 about BIA (SIA) effects dominating 
in large (narrow) gap systems must be understood as re- 
ferring to the dependence on the gap of the prefactors 7c 
and the part of an not proportional to the expectation 
value of the electric field. But-other factors such as the 
well width or an applied biasEj can change the relative 
contributions. 

Finally, the electron spins are also affected by the in- 
clusion of BIA terms. In Fig. |l^ the spins of the lowest 
conduction subband are shown for states lying on a circle 
in the kx — ky plane of radius k — 0.01 A. The direction 
of the spins has changed respect to the case without BIA 
terms (see Fig. ^l]). As it can be deduced from Eq. ([23[), 
it corresponds to the vector sum of the spins in Fig. |11| 
and the spins in Fig. each one weighted by their cor- 
responding splitting coefficients a's. 



V. CONCLUSION 

We have presented an implementation of an 8-band 
Finite Difference Effective Mass Approximation (EMA) 
method for calculating band structures. This implemen- 
tation is faithful to the Td microscopic symmetry of bulk 
zincblendes. We show that this allows the heterostruc- 
ture Hamiltonian to reproduce the D2d or the C21, point 
group symmetry of [001] grown zincblende heterostruc- 



tures. As a consequence, all symmetry effects close to 
the zone center, including the spin splitting of the sub- 
bands, are correctly described. A simple expression for 
the k^ splitting coefficient jc is given in terms of the 
k • p parameters. We find an extra contribution to 7c 
that can amount to a few percent of its value. When the 
method is applied to symmetric heterostructures, linear 
splittings in k are predicted as a consequence of the re- 
duced symmetry. This is not the case for standard EMA 
implementations. The bands of asymmetric heterostruc- 
tures are also studied and described in the context of the 
BIA model Hamiltonian for heterostructures. It is seen 
that, in the case studied, the SIA and BIA contributions 
to the spin splitting are of the same order of magnitude 
{an = 40.3 x 10-^° eV-cm and q;bia = 15.0 x 10-^° 
eV-cm respectively) even though the well is composed by 
narrow gap materials. Therefore, an accurate description 
of the bands will require the inclusion of both effects. 
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APPENDIX A: EXPLICIT FORM OF THE 
8-BAND k p HAMILTONIAN 

The 8-band k • p Hamiltonian derived by Trebin et 
al. ^ has been constructed using the theory of invari- 
ant£3, and thus it correctly describes the Td symmetry of 
bulk zincblendes. In particular, it includes terms break- 
ing the spin degeneracy of the bands at a general point in 
the Brillouin zone, making it ideal for the study of inver- 
sion asymmetry effects. It also accounts for the effects of 
strain and an external magnetic field. 

Although the Hamiltonian is explicitly shown in 
Ref. it is rewritten here for completeness and to cor- 
rect two typographical errors. In order to minimize the 
probability of introducing hard-to-detect errors in the 
definition of the matrix eleraciits, the Hamiltonian was 
first entered in Mathematical. From there, the appro- 
priate C code for each matrix element was generated au- 
tomatically with the instruction CForm. This method also 
has the advantage that it allows the algebraic operation 
of the Hamiltonian to find analytical forms for th e dis - 
persion relation very near the zone center (cf. Sec. |ll A| ). 

The 8-band k • p Hamiltonian will be expressed in the 

basis {|r6,+i), |r6,-i), |r8,+|), |r8,+i), |r8,-i), 

iFg,— |), |r7,-fi), jFy,— i)}. It can be written in a 
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H'"'=E, + E,+ '^+ A'e - gs^^a • H + Ci tre 

Hr=-^ {hifc' - 72 [(J.' - i J') fc^ + cp] - 273 [{J.J,} {k^ky} + cp]} - ^ { (kJ, + gj,3) + cp} 
Hr=Da trf + |D„ [(J^ - i J^) + cp] + |7:>; [2 { J, JJ + cp] 
Hm=^,C [{J. (J,^ - J|) } fc. + cp] 

JJ''==-A,o + ^jik" - 2k^^o- ■ H + i3d trf 

[P (fc,r, + cp) + iB (T, + cp) + iCa (r.cj,, + cp)] 

J/'="=-;7| [P (fc.Po; + cp) + iP (p;^ + cp) + iC2 (p^e^z + cp)] 

rrvs TTVS , ttvs i rjvs 

Hr=-^ [-372 + cp) - 673 (f/.y {k^ky} + cp)] - ^1 [U.H, + cp) 

Hr=2Du {U^^e^^ + cp) + 2D'^ {2U^y€^y + cp) 

Heh=^ [C4 (esy — tzz) fco; + C5 {txyky — t^zk^)] Ux + Cp 

cp means cyclic permutation, {AB} = + P^), tre — e^x + tj/y + tzz 

TABLE V: Matrix elements of the 8-band k ■ p Hamiltonian. 



O Td 

|r8,-3/2) |r8,+i/2) 

|r8,-i/2) |r8,+3/2) 

|r8,+i/2) ir8,-3/2) 

|r8,+3/2) |r8,-i/2) 



where, of course, (^H"^^ ~ H^", c refers to the two 
conduction band (CB) states, v to the four heavy hole 
(HH) and light hole (LH) states and s to the two spin- 
orbit split off (SO) states. 



TABLE VL Correspondence in basis functions for O and Td 
in Table 83 of Ref. 



block diagonal form 







JJCV 


H 




JJVC 




H 




\JJSC 


JJSV 


H 



(Al) 



The constituent blocks of the Hamiltonian are shown 
in Table The phases of the wavefunctions and the 
prefactors in Table ^ are chosen in a way that all the 
parameters are real. The cr and the p matrices are the 
Pauli matrices; the T matrices are given by 



rp _ 1 f -\/3 010 

^ 3\/2 I -1 V3 



T — 



-1 \/3 

3V2 V -V3 10 

-1 -\/30 

2V6 V V3 1 



VSO 1 Q 

3\/2 I 1 



rp _ 1 / -1 -V3 

3V2 1 \/3 1 



1 / -1 A 

2V6 V V3 -1 / 



T — V2 / 1 ^1 
— 3 (.0010^ 



zz 



3 vo -1 0; 

i / -1 ^ 
-1 , 



(A2) 



The matrices U are simply given by Ui —Tj. 

Note that a typographical error in the matrix T^x in 
Table I of the original article by Trebin et al£3 has been 
corrected here. The last equation in the group (A3) of 
Ref. ^ must also be corrected: 



X12 = -I (xf ^ - /2 (A3) 



APPENDIX B: BASIS STATE LABELS IN THE 
KDWS TABLES 



There is another remack about a point that can lead to 
confusion. Kostcr et a/.E3 have developed a set of tables 
(KDWS tables) for the Clebsch-Gordan coefhcicnts for 
point groups that are very helpful when constructing ex- 
plicit subspace-invariant matrices or when checking the 
symmetry properties of the Hamiltonian. However, in 
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their Table 83 the values they show can be used as dis- 
played for the O point group, but for T^i the values should 
be taken according to the lookup table for the basis state 
labels shown in Table [v| . 



APPENDIX C: FINITE DIFFERENCE METHOD 
SOLUTION OF THE EFFECTIVE MASS 
APPROXIMATION EQUATIONS 

In this appendix we describe the computational 
method used to find the hcterostructure bands and eigen- 
states. Follo»dng the Effective Mass Approximation 
(EMA) theoryB, the 8 band k-p Hamiltonian of Eq. ( [Ai| ) 
is transformed into a set of eight linear, second order, or- 
dinary differential equations. The appropriate boundary 
conditions are enforced and the equations are solved by 
means of a finite difference scheme. 



1. EMA Hamiltonian 



Starting from the usualElEaiia expansion of the wave- 
function ^k|| (r) as a linear combination of zone center 
Bloch states m„o modulated by a position dependent en- 
velope Fn{z) {n is the band index): 



M/k,, (r) = ^e"'«--F„(z)«„o (r) , 



(CI) 



the EMA prescription states that if ^fkn (r) is the solution 
of 

[H (k|| , k,;z) + U (z)] *k|| (r) = ^ ^k,, (r) , (C2) 

then Fn{z) will be the solution to 

[H (k|| , z)+U (z)] F(z) = E F(z), (C3) 

where k|| ~ k^x + kyij is the electron wavevector in the 
kx — ky plane, U (z) is an external potential that varies 



only in z (the change of material as a function of z is 
included in (ky, — zc^^; z)) and F{z) is a multicompo- 
ncnt vector constructed from the different _F'„(2:)'s. In 
an 8-band theory, F would have 8 components, each one 
multiplying the conduction band (CB), heavy hole (HH), 
light hole (LH) and split off (SO) basis states. 

On the other hand, the bulk k • p Hamiltonian can be 
expanded into its polynomial form for fc^ in the following 
manner: 



H (k) = H(2)A:2 + jid) (k||) + h(") (k|| 



(C4) 



Putting together Eq. ( |C3| ) and Eq. (C4), we must solve 
the following system of coupled differential equations to 
obtain the energies and eigenstates of the system: 



-H(2)52 - zH(i) (k||) + H(") (k||) + U{z)\ F{z) = 

EF{z). (C5) 
2. The Finite Difference Method 



There are several methods to solve numerically the sys- 
tem of c oupled ordinary differential equations |-giwn by 
Eq. (|C5|), such as the-transfer-matrix methodESlO, the 
finite element methodl23, the basis expansion methocO, 
etc. In this study, the finite difference method (EDM) 
has been employed because of its conceptual simplicity, 
the ease of introduction of arbitrary fields, its ability to. 
describe tunneling phenomena with only a few changesEEl 
and its numerical stability with respect to the transfer- 
matrix method, which requires the truncation of growing 
exponential statesEJ. 

In the finite difference method, the differential opera- 
tors are first written in a Hermitian form and then sub- 
stituted by their finite difference approximations over a 
discrete mesh (seCpE'lg. |5]) with N points. Following 
Chuang and ChangEJ, the following discretization scheme 
is used: 



J 



H(2)(z,+i) + H(2)(z,) 



i+l) 



2(Az)2 

H(^)(z,+i) + 2H(^)(z,) + H(2)(z,_i) 
2(Az)2 



, , H(2)(z._i) + H(2)(zO ^, 
f{zi) + 2{Az)^ (C6) 



-zH«(z)a,/ 



T 



H(i)(z)a,/-f a,HW(z)/ 

.H«(z,+i)+HW(zO 



4Az 

r 



f{z,+i)+l f[Z^-l) (C7) 



where Az is the separation between the mesh points, and Zi is the position of the i-th mesh point. 
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N ow, the application of the above equations to 



Eq. ( C5 ) yields the following system of N algebraic equa- 
tions: 



H, 



+ Hi^iFi + Hi_i+iFi-(_i — EF 



(C8) 



where is the eight- vector containing the envelope func- 
tion components corresponding to the i-th mesh point. 
This eigenproblem can be written in matrix form to bet- 
ter appreciate its sparse structure: 



Vh 



Ho,o 


Ho,l 











Ho,-i 


Hi,o 


Hii 

















H2,l 


H2,2 


H2,3 




















HAr_2,Af- 




HAr_2,Ar 


^N-l,N 

















EF 



(C9) 



where F is a column vector composed of the different 
F,'s. 

The resulting discretized 8x8 Hamiltonian matrices 
will be valid for both inner and interface mesh points 
without further modification: 



H, 



H 



(2) 
i+1 



2H 



(2) 



H 



(2) 
i-1 



Hf ' + [/, (CIO) 
(1) ^h(^) 



2(Az) 



H, 



H 



(2) 



H 



(2) 



H 



i+l 



H, 



H 



2{Az) 

(2) 
i-l 



4Az 



H 



(2) 



+ i- 



H 



(1) 



H 



(1) 



(Cll) 
(C12) 



2(Az)2 ' " AAz 
Hq.-i and Hjv-i,jv in Eq. (|C9| ) express the boundary 



conditions (BCs) of the problem. When studying a quan- 
tum well, the BCs arc that the wavefunction must vanish 
far from the well region. This is accomplished by setting 
the barrier region wide enough, and requesting 



r_i 



•n 



which translates into 



(C13) 



N1/N2/N3 Quantum Well or Superlattice 



Ho. 1 H 



N-1,N — 



(C14) 



Monolayer! 
Boundary < 
Points I 



Mesh 
Points 



1 





N,-l N, 



N,-l N, 



When finding the energies and states of a superlattice, 
the Bloch BCs apply, and the envelope function is re- 
quested to have the supercell periodicity d, modulated 
by a phase: 



Fn = e''«''Fo 



F_i = e^'^^F 



Hat-i^At — e*^''Hjv_2,7V-l 

(C15) 
(C16) 



FIG. 15: Schematic of the structure under study, with points 
separating the monolayers, and mesh used when solving the 
Effective Mass Approximation equations. The mesh points 
need not coincide with the monolayer boundaries. 



where q is the electron wavevector along the z direction, 
and it has been assumed that the same material is at 
mesh points and iV — 1 . 
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